An all-optical event horizon in an optical analogue of a Laval nozzle 
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Exploiting the fact that light propagation in defocusing nonlinear media can mimic the transonic 
flow of an equivalent fluid, we demonstrate experimentally the formation of an all-optical event 
horizon in a waveguide structure akin to a hydrodynamic Laval nozzle. The analogue event horizon, 
which forms at the nozzle throat is suggested as a novel platform for analogous gravity experiments. 
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Event horizons are well known in the context of as- 
trophysics and cosmology. Less known are the analogy 
between an astrophysical event horizon and the sonic 
horizon in transonic fluid flow, and the prediction that 
a thermal spectrum of sound waves should be emitted 
from a sonic horizon, in analogy with Hawking radia- 
tion [1]. These analogies set the stage for attempts to 
create laboratory black hole analogues, involving various 
physical scenarios, from water flowing in a channel to the 
acceleration of a superfluid to nonlinear optical experi- 
ments [2l4Io|. We have recently proposed an alternative 
approach to analogue gravity experiments - an all-optical 
experiment based on laser light propagation in a distinc- 
tive nonlinear waveguide, which is analogous to a Laval 
nozzle (a well-known device in the context of aerodynam- 
ics). This approach has two great advantages over previ- 
ous experiments: The attainment of supersonic velocities 
is very easy, and the analogue of Hawking radiation has 
a unique optical signature, which can be readily detected 
The analogy is based on the realization that under 
certain conditions light can "flow" through certain types 
of media in a fashion reminiscent of actual fluid flow. A 
prime example is a laser beam propagating through a 
Kerr-type nonlinear medium, which is usually described 
analytically by the Non-Linear Schrodinger equation 
The latter can be mapped, through the Madelung trans- 
formation [12[, to a pair of coupled equations for the 
amplitude and phase, which have the form of contin uity 
and Euler equations for an equivalent fluid 0, 0, 
which may be called " luminous fluid" : 



d z p + V- [pv] =0, 
d z v + \vv 2 = -tt V 
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Eqs. (pQ) describe the evolution of the complex ampli- 
tude f(x,y)e~ t(p ( x,y " ) as the light propagates along the z 
axis with the wave vector /3q. Here p = f 2 is the light in- 
tensity, and the transverse component of the wave- vector, 
Pov = — V<£, plays the role of velocity. The coordinate 
z plays the role of time, /3q is equivalent to the mass 
of a particle, and the spatially inhomogeneous refraction 
index assumes the role of a potential, U ex t(x,y). The 



FIG. 1: (Color online) Images of the waveguide structures, (a) 
shows the input plane, with six circular openings of lightpipes 
of different dimensions, (b) shows the exit plane, with grooves 
of triangular cross-section, cut half-way along the sides of the 
channels, and forming the divergent sections of the nozzles. 



nonlinear term is due to the Kerr effect. Thus incident 
light, which propagates at an angle relative to the z axis 
is mapped onto a fluid with a finite transverse velocity, 
and a change of that angle corresponds to acceleration of 
the fluid. This approach has proved to be an extremely 
powerful one when applied to the problem of coherent 
tunneling [l8l-[23|. It has also been used to model dis- 
persive shock waves that appear when the nonlinearity is 
repulsive (i.e. self-defocusing, A > 0), and consequently 
an equivalent real sound velocity s 2 = Xp/Po can be de- 
fined [T3-[l7l, |24] - [26| . Our proposal for an optical ana- 
logue of the Laval nozzle uses a similar approach, [gl [2?J 
and shows that when the light is confined by a prop- 
erly shaped channel (i.e. waveguide), represented by the 
transverse potential U ex t > it may propagate transversally 
in a way that resembles the accelerating transonic flow 
of the equivalent, luminous fluid. Thus low-incidence, 
(i.e. "subsonic") laser light is predicted to accelerate, 
(i.e. change its propagation direction) while traversing 
the nozzle, reaching a critical velocity, which is equiva- 
lent to the sound velocity in a real fluid, at the nozzle 
throat, and exiting the nozzle at a "supersonic" velocity. 
While the notion that light propagation in nonlinear me- 
dia can mimic the flow of an equivalent fluid is not new, 
the equivalent of transonic acceleration in an optical ana- 
logue of a Laval nozzle has never been observed before in 
an actual experiment, and is demonstrated here for the 
first time. 
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FIG. 2: (Color online) Waveguide exit plane images and cor- 
responding power density cross-sections along the nozzle axis 
for an input power of 2 Watts and an iodine concentration of 
~40 ppm. (a) and (b) show data for a 2 mm diameter waveg- 
uide, (c) and (d) show data for a 3 mm diameter waveguide, 
and (e) shows the free expansion of the beam outside the 
waveguide structure, (f) is a reference image of the 3 mm 
diameter waveguide. 



We study the flow of luminous liquid through an opti- 
cal Laval nozzle experimentally by launching a continu- 
ous wave laser beam into an appropriately shaped waveg- 
uide with reflective walls, filled with a Kerr- type defocus- 
ing nonlinear material. The experimental challenge here 
is to create conditions of steady flow with a subsonic in- 
put velocity. Since the velocity is k x //3o (where k x is the 
x-component of the wave vector) and the sound velocity 
is s 2 = \p/Po, the above input conditions imply a small 
input angle of the beam and a high nonlinearity and/or 
input intensity (Note that, for a given angle, low and high 
intensities correspond to supersonic and subsonic flow, 
respectively). However, an unavoidable consequence of 
these conditions is strong self-defocusing of the beam, 
and as a result the wave packet, which traverses the noz- 
zle is an expanding "droplet" of liquid, with a tendency 
of the power density in the cavity to decrease with in- 
creasing input power. Furthermore, while the peak in- 
tensity of the droplet may correspond to subsonic flow, 
it is always surrounded by supersonic flow (in contrast 
to the usual case in hydrodynamics), and when confined 
to a Laval nozzle such as the one discussed in Ref. 0, 



the fluid flows from the throat towards both sides of the 
nozzle. It is thus impossible to generate the steady sonic 
background flow conditions stipulated by the theory in 
a simple waveguide with a convergent-divergent cross- 
section formed by two convex walls. To circumvent this 
problem we use an alternative waveguide design, based 
on a light-pipe of circular cross-section drilled in an alu- 
minum block, with a groove of triangular cross-section, 
cut along the side of the channel, acting as the divergent 
section of the nozzle (see Fig. [I]). The total length of 
each light pipe is L = 67 mm, and the groove extends 
over the second half of this length. This design is in- 
tended to "trap" the expanding beam and confine it in a 
homogeneous, high-density and low- velocity mode, thus 
preparing it for ejection through the groove, and is akin 
to the configuration of a rocket engine: a high-pressure 
gas is first loaded into a combustion chamber, and is 
then expelled through the nozzle. The aluminum block, 
with several nozzles of different diameters, aperture sizes 
and groove opening angles, is enclosed in a plexiglass 
cell with glass windows, which is sealed and filled with 
iodine-doped ethanol. The nonlinear index variations re- 
sult from optical absorption by the iodine, which in turn 
leads to thermally-induced changes of the index of re- 
fraction - a non-local nonlinearity, which slightly washes 
out the thermal gradients |17j . The nonlinearity Xp can 
be expressed, in terms of the nonlinearly-induced refrac- 
tive index change 5n, as 5nf3o/no, where no is the linear 
refractive index of the material [17j. The corresponding 
dimensionless sound velocity is then s 2 = Sn/no : mean- 
ing that the input beam is subsonic for k x //3o < ^Sn/no. 
We use a continuous-wave frequency-doubled YAG laser 
(532 nm), and focus the beam to a ^ 0.5 mm waist at 
the input of a waveguide. The input power is varied by 
means of the laser controller, in order to avoid thermal 
effects in variable-density filters. Images of the exit plane 
of the waveguide are recorded by means of a CCD cam- 
era. In all cases images were acquired after stabilization 
of the thermal gradients. 

Figure [2] presents images of the exit plane of two of 
the waveguides, and the corresponding power density 
cross-sections along the nozzle axis, for an input power 
(2 Watts) that is sufficiently high to completely fill the 
waveguides (at an iodine concentration of ~ 40 ppm). 
Figure EJa) shows a 2 mm diameter waveguide, and Fig. 
[2](c) shows a 3 mm diameter waveguide, both having a 
~ 0.5 mm opening (i.e. nozzle throat). Figures [2fb) and 
[2]^d) are the power density cross-sections corresponding 
to Figs. [2fa) and[2fc), respectively, obtained by sum- 
mation over 12 CCD lines at the center of each noz- 
zle. Figure [2^e) shows the free expansion of the beam 
when it propagates outside the waveguide structure, and 
Fig. [2ff) shows a reference image of the 3 mm diame- 
ter waveguide, obtained with incoherent light and with 
the laser beam blocked. Figures [2fa)-(d) clearly show 
the jets of luminous liquid ejected from the nozzles as 
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the beam propagates through the waveguides. Note that 
the jets extend farther than the edge of the beam un- 
dergoing free expansion (Fig. [2fe) - A detailed analysis 
is presented in Fig. [3j). Furthermore, there is a sharp 
drop in the density as the jet exits the nozzles, which is 
clearly seen in the images and in the power density cross- 
sections. This demonstrates that the luminous liquid is 
accelerating at the nozzle throat rather than gradually 
expanding through the opening. Finally, while the con- 
fined beam propagates along the waveguide walls at a 
very slow (i.e. subsonic) velocity, the following analysis 
shows that the jet of luminous liquid is indeed supersonic: 
The dimensionless velocity of the jet outside the waveg- 
uide is first calculated from its extension in the trans- 
verse direction, deduced from the images. The relation 
is simply v = Ax/Az = 2Ax/L ~ 0.1, where Ax is the 
transverse distance from the nozzle throat to the edge of 
the jet, and Az = L/2 is the distance along the z axis 
that the same part of the jet has propagated by the time 
it reached the exit plane. This velocity should be com- 
pared to the local sound velocity, which can be estimated 
by analyzing the light intensity distribution in the exit 
plane and the rate of expansion of the freely-expanding 
(i.e. self-defocusing) beam. The latter, deduced from 
Fig. Efe), allows us to calculate Xpo and the correspond- 
ing sound velocity at the input. The former in turn allows 
us to deduce the sound velocity, which corresponds to the 
lower density of the jet, taking into account the expan- 
sion of the beam in the light-pipe, the relative intensities 
of the jet and inside the light-pipe, and measured losses. 
This calculation gives a local sound velocity in the jet on 
the order of 1 x 10 -3 or less, meaning that the local Mach 
number is >100. This clearly establishes that the lumi- 
nous liquid undergoes transonic acceleration and forms a 
"sonic" horizon as it expands through the nozzle. 

In Fig. [3] the dimensionless velocities (v — 2Ax/L) of 
the jets emanating from the 3 mm and 2 mm nozzles are 
plotted as a function of the input intensity, for a fixed io- 
dine concentration. Also shown is the velocity at the en- 
velope of the freely-expanding beam, which we estimate 
as dx/dz « Ax' / L (In this case we measure Ax' from 
the center of the beam, which we determine from low- 
intensity measurements; This velocity corresponds to the 
asymptotic expansion angle, obtained for L » l/Apcb 
i.e. when the propagation distance is much longer than 
the characteristic defocusing distance, and is a reason- 
able estimate for intensities > 1 Watt). Figure [3] clearly 
shows that the velocity of the jet is higher than that of 
the freely-expanding beam throughout the experimental 
intensity range, in spite of the fact that the initial condi- 
tions for the free expansion involve higher pressures (the 
equivalent of pressure in a luminous liquid is P = ^\p 2 ). 
On the other hand Fig. [3] shows that the jet emanating 
from the 2 mm waveguide is slower than that ejected from 
the 3 mm nozzle, although, for a given input intensity, 
the pressure in the latter is supposed to be lower. This 
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FIG. 3: The measured jet velocities as a function of input 
power for an iodine concentration of ~ 40 ppm. Data are 
shown for two waveguides and for free expansion of the beam 
outside the waveguide structure. The calculation of velocities 
is explained in the text. The lines are guides to the eye. 



discrepancy can result from the fact that the opening in 
the 2 mm waveguide subjects a larger angle, resulting 
in a less directional jet (compare Figs. [2ja) and[2fc)). 
We also measure higher losses (due to scattering and ab- 
sorption) in the smaller waveguide, so in fact the power 
densities in the two waveguides are comparable. Finally, 
the non-locality of the nonlinearity may have a stronger 
effect on the 2 mm waveguide. 

Measurements at lower intensities illustrate another 
regime of operation of the nozzles. Figure ffl[a) shows 
an image of the 3 mm waveguide for self-defocusing (i.e. 
beam expansion) that is not sufficiently strong to com- 
pletely fill it (an input power of 2 Watts and an iodine 
concentration of only ^20 ppm). Figure Hfb) shows the 
corresponding power density cross-sections along the noz- 
zle axis. In order to produce the jet seen in Fig. SJa) the 
nozzle throat had to be displaced (horizontally) relative 
to the beam axis. Figure |4^c) shows the dependence of 
the jet velocity on the displacement and the input power. 
Note that as the input power increases the optimum ac- 
celeration is obtained when the beam axis is moved far- 
ther away from the throat (at an input power of 2 Watts 
the beam axis is then near the center of the waveguide). 
A comparison of the data in H^c) with the divergences 
of the freely-expanding beam, measured separately for 
the same input powers, shows that optimum accelera- 
tion at the nozzle is obtained when the envelope of the 
freely-expanding beam coincides with the nozzle throat 
halfway through the waveguide (i.e. at z = L/2). Under 
these conditions the luminous liquid entering the nozzle 
is already supersonic, and is accelerated further in the di- 
vergent section of the nozzle. The smooth power density 
cross-section shown in Fig. SJb) supports this interpreta- 
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FIG. 4: (Color online) Acceleration of supersonic flow, (a) 
An image of the exit plane of the 3 mm waveguide, (b) The 
corresponding power density cross-section along the nozzle 
axis, (c) The jet velocity as a function of the distance between 
the beam axis and the nozzle throat, for four input powers 
(The curves are guides to the eye). The iodine concentration 
is ~ 20 ppm. 



the throat. This ratio needs to be of order unity so that 
both parts would be observed in the experiment, and 
allow direct measurement of the Hawking temperature 
Tjj- The same condition can also be written as Ih « 2/q, 
where Ih = fi/Tn = ^cj\fZ~s and Iq = l/v$ are charac- 
teristic length scales of the horizon and the fluctuations, 
respectively. In the experiment described here Ih is on 
the order of a few meters (c~ lx 10 -3 m, s ~ 1 x 10 -3 ), 
while lo must be on the order of a few centimeters (the 
length of the cavity). This may still allow observation of 
the part of a fluctuation which is carried away with the 
supersonic flow, but the part which penetrates into the 
subsonic region will most likely be submerged in noise. 
Therefore Ih must be decreased by two orders of magni- 
tude. Note, however, that for given input intensity and 
nonlinear coefficient the factor c/s in the expression for 
Ih grows like c 2 . We therefore estimate that an order of 
magnitude smaller cavity will be sufficient for observing 
both parts of a straddled classical fluctuation. The re- 
quirement for a slow rate of acceleration can be met by a 
refined, smoother waveguide cross-section, compared to 
the rudimentary prototype that we have used here for 
demonstration purposes. 

We thank Mr. G. Elazar for his technical support. 
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tion (compare this to the sharp density gradients at the 
throat in Fig. [2j). In this case the nozzle operates in a 
regime that is not typical of a Laval nozzle, though. 

In conclusion, we demonstrate experimentally, for the 
first time, the transonic acceleration of a luminous liq- 
uid through an optical analogue of a Laval nozzle. The 
analogue of a sonic event horizon, which forms at the 
nozzle throat, lends itself to studies of classical (and pos- 
sibly quantum) fluctuations that are akin to Hawking 
radiation Q. Compared to other experiments and pro- 
posals for "analogue gravity" llOL our experiment 
has the advantage that it allows easy generation of su- 
personic flow conditions. While the nonlinearity length 
l n i = 1/VApoA) ~ 30/im in our experiment is sufficiently 
short to give way to fluctuations with a linear disper- 
sion relation, [9] the challenge is to create an equivalent 
Hawking temperature that is high enough to measure 
experimentally. Note that this is not a real tempera- 
ture, but it rather a constraint on the minimum light 
intensity (and sound velocity) required in the waveguide: 
As explained in Ref. 9, the ratio of amplitudes of the 
two parts of a classical fluctuation with wave vector v Q - 
the part which is carried away with the supersonic flow 
(i.e. "falls" into the black hole) and the part which pen- 
etrates into the subsonic region (i.e. "escapes" from the 
black hole), is exp(^|^) > 1, where c is a characteristic 
length scale of the nozzle and s is the sound velocity at 
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